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^ ' Pionium lifetime corrections to the nonrelativistic formula are calculated in the frame- 

^ I work of the quasipotential-constraint theory approach. The calculation extends an earlier 

I evaluation, made in the scheme of standard chiral perturbation theory, to the scheme of 

O ■ generalized chiral perturbation theory, in which the quark condensate is left as a free 

^ ! parameter. The pionium lifetime is calculated as a function of the combination (oq — Oq) 

r-| ■ of the UTT 5- wave scattering lengths with isospin / = 0, 2. 

PACS numbers: ll.30.Rd, 12.39.Fe, 13.40.Ks, 11. 10. St. 

Oh! 



Keywords: Chiral perturbation theory, Pion, Hadronic atoms. Electromagnetic cor- 



^ ' rections, Relativistic bound state equations. 



The DIRAC experiment at CERN |jl| is expected to allow, in the near future, the 
measurement of the lifetime of the pionium (tt+tt" atom) with a 10% precision. The latter, 
in turn, through its relationship with tttt scattering lengths, would provide a determination 
of the combination (oq — Oq) of the S-wave scattering lengths with isospin / = 0, 2 with 
5% accuracy. The strong interaction scattering lengths Oq and have been evaluated 
in the literature in the framework of chiral perturbation theory (xPT) to two-loop order 
of the chiral effective lagrangian P, ^, Q . Therefore, the pionium lifetime measurement 
provides a high precision experimental test of chiral perturbation theory predictions. 

The nonrelativistic formula of the pionium lifetime in lowest order of electromagnetic 
interactions was first evaluated by Deser et al. ^ and later reanalyzed by others 0. It 
reads: 
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where ip^ (0) is the wave function of the pionium at the origin (in x-space). 

An evaluation of the relativistic and higher-order electromagnetic corrections to the 
above formula was recently done by several authors. In the frameworks of quantum field 
theory and xPT, three different methods of evaluation have led to similar estimates, of 
the order of 6%, of these corrections [0, H, The first method uses a three-dimensionally 
reduced form of the Bethe-Salpeter equation (the quasipotential-constraint theory ap- 
proach) and deals with an off-mass shell formalism |^. The second method uses the 
Bethe-Salpeter equation with the Coulomb gauge 0. The third one uses the approach 
of nonrelativistic effective theory [y, |TU|, [IT], |TB[. (A survey of other approaches is 
presented in the second paper of Ref. 0.) 

The pionium lifetime, with the sizable relativistic and electromagnetic corrections 
included in, can be represented as: 

i = r = 

r 
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where 7?.e7Vloo,+- is the real part of the on-mass shell scattering amplitude of the process 
7r~^TT~ 7r°7r°, calculated at threshold, in the presence of electromagnetic interactions 
and from which singularities of the infra-red photons have been appropriately subtracted 
14| ; the factor 7 represents the contribution of interactions at second-order of perturba- 




tion theory with respect to the nonrelativistic zeroth-order Coulomb hamiltonian. The 
explicit expressions of 7?.eA^oo,+- and of 7 may differ from one approach to the other, 
according to the way the singularities of the infra-red photons are subtracted, but their 
total contribution should be the same. 

The purpose of the present article is to calculate the pionium lifetime in the framework 
of generalized xPT []T3|, ^. In this framework, according to the observation that the 



fundamental order parameter of spontaneous chiral symmetry breaking in QCD is F-j^, the 
decay coupling constant of the pion, the quark condensate in the chiral limit < 0|gg|0 >o 
is left as a free parameter. Its value should depend on the details of the mechanism of 
chiral symmetry breaking and would require an independent experimental test. In the 
framework of standard xPT, the mechanism of chiral symmetry breaking would be very 
similar to that realized in a ferromagnetic medium [0. In this case, the value of the 



Gell-Mann-Oakes-Renner (GOR) parameter [T?!, defined as 



2m < 0|gg|0 >o 
^GOR ' 

TT TT 

where 2rh = niu + trid and is the pion physical mass, would be close to one. Equiva- 
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lently, the quark condensate parameter 

^^_<0||0>,_ 

where F is in the chiral SU{2) x SU{2) hmit, would be of the order of the hadronic 
mass scale ~ 1 GeV. This assumption fixes the way standard xPT is expanded: 
the quark condensate parameter B is assigned dimension zero in the infra-red external 
momenta of the Goldstone bosons, while the quark masses are assigned dimension two 

i- 

On the other hand, in an antiferromagnetic medium, one meets a situation where the 
quark condensate would be zero or very small |I6| , |T8[. Recently other possibilities were 



also advocated ||T9|, : the existence of a possible chiral phase transition in QCD [^T], |22 



at relatively low values of the light quark flavor number might induce, in the standard 
case, a strong flavor dependence of the quark condensate, which would have the tendancy 
to decrease by passing from SU{2) x SU{2) to SU{3) x SU{3). 

The framework of generalized xPT offers the possibility of experimentally testing 
various issues of chiral symmetry breaking mechanism. In this framework one relaxes the 
GOR assumption and treats the order of magnitude of the quark condensate parameter 
i? as an a priori unknown quantity (awaiting a precise experimental information about 
it) leaving to it the possibility of reaching small or vanishing values. To this aim, B is 
assigned dimension one in the infra-red momenta of the external Goldstone bosons and 
accordingly quark masses are also assigned dimension one. Due to this rule, at each order 
of the perturbative expansion, generalized xPT contains more terms than standard xPT. 
For instance, the pion mass formula becomes at leading order: 

(m^)g = 2mB + 4mM, (5) 

where the constant A is expressible in terms of two-point functions of scalar and pseu- 
doscalar quark densities. In the standard xPT case, this term is relegated to the next- 
to-leading order. 

Processes involving only Goldstone bosons are sensitive at leading order to the value 
of the quark condensate parameter B. Among them, the vrvr scattering amplitude plays 
a key role. At the tree level of the chiral effective lagrangian, the amplitude A{s\t,u) has 
the expression: 

A{s\t,u) = ^{s-2mB), (6) 

which displays explicit dependence on the quark condensate parameter. 

At the one-loop order the strong interaction on-mass shell scattering amplitude can 
be described by four parameters, a, (3, Ai and A2, and takes the following form p3| : 

A{s%u) = M',l- = 4-Js - -ml) + a ^ 
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where m^^ is the physical pion mass, J the conventional loop function [0] and Am? A the 
quadratic mass term present at the tree level in ml [Eq. (^]. The parameters Ai and 
A2 are related to the standard renormalized low energy constants l\ and I2 0; their 
expressions (at the one-loop order) are: 

Ai = 211 
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(/i is the renormalization mass.) j3 is essentially related to the deviation of from F. a 
is mainly related to the quark condensate parameter; its expression at the tree level is: 

Am'A 
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The one- loop expressions of F', ml, (3 and a are given by the following relations: 
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where are combinations of renormalized low energy constants present in the generahzed 
xPT lagrangian 0, |2^, and l"^ and II are the renormahzed low energy constants already 
present in the standard x^T lagrangian 0. [Our effective lagrangian slightly differs from 
that of Refs. |4| in that we continue using in it the /4-term, instead of replacing it 
through the equations of motion with the ^^^-'-term.] 

In terms of the above parameters the expressions of the S- and P-wave scattering 
lengths are: 

^ \/3-^!gl(Ai-A2)+ , , ^(5a^-40a/3-16/j^). (17) 



At the two-loop order, the vrvr scattering amplitude is described by six parameters, 
a, /3, Ai, A2, A3 and A4. The five parameters other than a are weakly dependent on the 
quark condensate; thus a remains the only parameter to be strongly sensitive to the quark 
condensate value. It is evident from the tree level relation (0) that in the standard case 
a remains close to 1, while in the extreme case of generalized xPP (antiferromagnetic 
case), where the quark condensate vanishes {2rhB ~ 0), it approaches the value 4. The 
expressions of the scattering lengths and of the effective ranges at the two-loop order in 
terms of the above six parameters can be found in Ref. 0]. 

In order to obtain the pionium lifetime, one first needs to calculate 7?.eA1oo,+- [Eq. 
(^)]. This includes, in addition to the strong interaction contributions, those of the 
electromagnetic interactions. These effects are taken into account in the chiral effective 
lagrangian by the presence of new terms involving low energy constants p5| , [26| , In 
generalized xPT the effective lagrangian contains more terms than the standard one, these 
being proportional to the low energy parameter A [Eqs. (P) and (|10])]. The additional 
elctromagnetic terms needed for the present process are given, in standard notations, by 
the following part of the effective lagrangian (in the SU{2) x SU{2) case): 

Z:^^^^^) = \F'A[h,{QUQU^){x'U + U^x? 

+8he{QUQU^){x^x) + h7{Q'){x^U + U^xf 

+hMQUx^U + U^xU^qY) + ho{QUx^U){U^xU^Q)], (18) 



where Q is the quark charge matrix, Q = exdiag.(2/3, —1/3), and the sources x have been 
restricted to a combination of an isosinglet scalar density and an isotriplet of pseudoscalar 
densities, x = + ^P 't, with s and p real, r representing the Pauli matrices. The renor- 
malization of the coefficients k is done, with dimensional regularization in d-dimensional 
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space-time, according to the decomposition 

h = + Kif^), A = - ^(ln(47r) + r'(l) + 1)), (19) 

fj, being the renormahzation mass. The /^-dependence of the renormahzed coefficients 
kl{n) is fixed by the prescription that the ks are /^-independent. The coefficients k are 
found to have the following values: K15 = | + 6Z, kiq = —Z, Kn = — ^Z, K19 = 2, 
K20 = 5, where Z = C /F^, C being the coefficient of the lowest-order, 0{e^p^), electro- 
magnetic term 

^(eV) = dQUQU^), (20) 
responsible for the pion mass shift at that order [p^ : 

{Aml)o = 2e'^, Ami ^ ml^ - mlo . (21) 

We use in the following for the other coefficients k the notations of Ref. . 

Treating the electromagnetic interaction in 7^67^00,+- as a perturbation yields for 
the zeroth-order part of the latter the strong interaction amplitude, which is already 
calculated in generalized xPT to two- loop order ^ . It is therefore sufficient to calculate, 
up to one-loop order of the chiral effective lagrangian, the first-order electromagnetic 
correction to it. It should be emphasized here that, at the numerical level, the strong 
interaction amplitude is calculated in the literature 0, 0, Q by identifying the pion mass 
appearing in it with the physical mass of the charged pion. 

The electromagnetic corrections can be divided into two categories. The first one con- 
cerns corrections that arise from conventional pion-photon interactions. The second one 
concerns corrections that arise from quark-photon interactions which manifest themselves 
in the chiral effective lagrangian through the presence of the 0(e^p°) mass shift term of 
the charged pion [Eq. (^)], which survives the chiral limit and produces the main part 



of the pion mass difference |27|. (It also induces, through renormahzation, higher-order 
counterterms in the chiral effective lagrangian.) 

The pion-photon interaction yields in general an infra-red divergent amplitude on the 
mass-shell. Usually, this divergence is avoided by giving the photon a small mass and 
then subtracting the singular pieces |T^. In our approach, based on the quasipotential- 



constraint theory method we use an off-maas shell formalism, where the total energy 
of the two-pion system is fixed at the bound state energy, that is, below thereshold. In this 
case the photon field can be taken massless and the corresponding scattering amplitude is 
finite. However, it still contains a spurious finite infra-red term (with respect to the bound 
state problem). This term is cancelled by the presence of three-dimensional diagrams, 
called constarint diagrams, the role of which is to trasnform, through the Lippmann- 
Schwinger equation, the scattering amplitude into an irreducible kernel or a potential. 
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The sum of the two contributions is finite and free of spurious terms |^. It then can be 
continued, neglecting O(e^) terms, to the threshold of the vr+vr" system. 

The quark-photon interaction terms are free of infra-red singularities and the corre- 
sponding part of the scattering amplitude can be readily calculated, taking into account 
in particular the pion mass difference. 

The electromagnetic corrections to the neutral and charged pion masses are given to 
order e^p^ by the following formulas: 

TV 
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(23) 



where 771^^ is the strong interaction mass (|T^. 

We present in the following the difference, A7^eA^oo,+-! at the vr+vr" threshold, of 
7?.eA1oo,+- [Eq. (H)] calculated in the presence of electromagnetism including the terms 
of order e^p° and e^p^ and of the strong interaction amplitude TZeM.QQ^_ [Eq. (^)] 
calculated with the physical charged pion mass: 
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+ ^(^)'[5(A;; + k{,) - 5{kl, + kl, + k[,) + 2k;,]}. (24) 

Here (a)o represents a in the chiral SU (2) x SU (2) limit [Eq. ([T0|)]; and correspond 
to the strong interaction quantities [Eqs. ([TT| ) - ([12|) ]; and are low energy constants 
appearing in the generalized xPT lag rangian p4|| ; At/i^ = vp? — ttl^o. 

We list below the remaining corrections to the pionium decay width |^. (We desig- 
nate by a^"^ the fine structure constant in order to distinguish it from the parameter a 
introduced previously.) 

1) The strong interaction correction coming from second order perturbation theory 
with respect to the bound state wave equation is: 

(Ar),,,, = 1.5a^"^(2a° + a^)ro. (25) 

2) The vacuum polarization correction is: 

(Ar)_. = 0.41a^-ro. (26) 

3) An effective 0{e^p^) correction also arises from formal O(e^p^) effects, corresponding 
to diagrams with one pion loop and one photon propagator. In the present bound state 
formalism, the diagram with one pion loop with one photon exchange provides an infra-red 
logarithmic contribution, which is: 

(Ar)o(,v) = -^(2a°-fa2)(21na'^- + 31n2 + 21C(3)/(27r2))ro 

= 2.2«^'"(2a° + a2)ro. (27) 

This effect could be considered as being part of 7^eA^oo,+- [Eq. (0)]; however, since 
A7^e7Vloo,+- in Eq. ( P^ has been defined as including only one-loop graphs, we write it 
here separately. 

4) The corrections due to isospin breaking in the quark masses being quadratic in 
(m„ — rrid) , provide negligible effects and are ignored. 

5) The kinematic correction coming from the phase space factor [Eq. (|^)] is not 
included in the definition of AF/Fo and should be directly incorporated in F or r. 

Collecting all these contributions, one finds for the total dynamical correction to the 
nonrelativistic decay width formula [Eqs. (^-(0)] the following expression: 

AT ^ 6A^^AW + 3.7^e.(2a° + ag) + 0.41a^-, (28) 

where af] and a, are the strong interaction scattering lengths calculated up to two loops 
in the framework of generalized xPT |Q. 
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In order to evaluate the pionium lifetime, one first needs to know the value of the 
combination (oq — Og) of the scattering lengths, which enters in the nonrelativistic formula 
(|I|) as well as in the correction (^). The evaluation of the latter necessitates also the 
knowledge of the other parameters of the theory. Among the six parameters [a, (3, Aj, 
i = 1, ... ,4) of the strong interaction amplitude five of them {f3 and the As) are weakly 
dependent on the quark condensate value. It is then natural to fix them at some mean 
values and to consider a as the only variable in the preceding relations. The precise values 
of the six parameters depend on the behavior of the tt — tt scattering amplitude in the 
whole low energy kinematic region. A determination of these values in the standard xPT 
case was presented in Ref. [^. In the present problem, however, the main quantity of 
interest is the combination (oq — Oq) of the scattering lengths, rather than the general 
scattering amplitude itself. Our aim is to establish a relationship between the pionium 
lifetime and this combination of the scattering lengths which could serve us to extract 
the value of (oq — a^), within an uncertainty interval, from the experimental value of the 
pionium lifetime. Furthermore, we want to place the prediction of standard xPT on the 
central line of the above relationship. 

To achieve this, we have fixed the values of the parameters (3 and Aj = 1, . . . , 4) from 
the predictions of standard xPT 0. We have considered set I of threshold parameters, 
for which one has in particular (ag — ag) = 0.258, and determined the parameters from the 
values of Oq, Oq, a\, a^, feg ^i- We have obtained the following values for the parameters: 
a = 1.021, P = 1.109, Ai = -10.28 x lO^^, Aa = 15.90 x lO^^, A3 = 0.81 x 10-^ 
A4 = —1.00 X 10~^. These values should be considered as a convenient means of analyzing 
the present problem, rather than the best values for the entire scattering amplitude; their 
role is to introduce the standard xPT prediction as an initial data placed on the central 
line of the relationship between the lifetime and (ag- ag)- Once the five parameters (3 and 
A are fixed, a one-to-one correspondence is established between a and the combination 
(oq — Oq) of the scattering lengths. Instead of considering a as the main variable of the 
problem, it is preferable to consider ( the principal variable of interest, since it 

has a direct physical meaning. We have thus considered a variation of (af] — Oq) within the 
interval 0.250-0.370 which essentially covers the domain of generalized xPT, including 
standard xPT. For each given value of (oq — Oq), the other combination (2aQ + Oq) is 
calculable through the low energy threshold parameter formulas [Q. 

The correction ( p8|) to the decay width depends also on the electromagnetic low energy 
constants k^. Some of these have been calculated in the SU{3) x SU{3) case in standard 
xPT in Ref. ||28[ and could be converted to the SU (2) x SU (2) case. We have assumed 
that the combination of the k^s present in the standard xPT case (the terms in Eq. ( p^ 
not proportional to (a — /3) or to {a — jS)"^) keeps its numerical value also in generalized 
xPT. As to the additional low energy constants (strong and electromagnetic) appearing 
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in generalized xPT, we have considered them as uncertainties of the order of l/(167r^) 
and added them quadratically. Furthermore, numerically (a)o in Eq. ( pSf ) has been taken 
equal to a. Finally, the following numerical values have been used: = 93.2 MeV, 
= 139.57 MeV, m^o = 134.97 MeV. We have taken m.^, which appears in higher 
order corrective terms, equal to and yU = rrip = 770 MeV. [An error exists in the 
conversion formulas of Ref. [0 of the combination of k^s of the standard x^T case, due 
to missing factors coming from the strong interaction low energy constants U and from 
the difference of chiral limits of in the SU{3) x SU{3) and SU{2) x SU{2) cases. 
The correct formula has been given in Ref. 0. Numerical predictions already made are 
however almost insensitive to this modification.] 

The uncertainty in the relative correction (|28|) was estimated in Ref. in the stan- 
dard xPT case, to be of the order of 2% with respect to Fq. This uncertainty has been 
maintained here and the additional one coming from the generalized xPT coefficients 
added to it. We have found that the uncertainty varies from 2% (standard xPT case) to 
2.5% (extreme case of generalized xPT). 
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Table 1: The pionium lifetime r as a function of Oq — Og, with the low energy constants (3 
and A fixed from the set I solution of standard xPT 0. tq is the lifetime obtained from 
the nonrelativistic formula. 

We present in Table |I] values of the pionium lifetime for some typical values of (oq — Oq), 
as well as the corresponding values of and the relative correction corresponding 

to Eq. (|8|). 
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We observe that the relative correction |Ar/ro| decreases with increasing a. This is 
mainly due to the fact that the tree level correction tends to zero as a approaches the 
value 4. 

In Fig. we have represented the curve of the lifetime r as a function of the combina- 
Oq) of the scattering lengths (full line). The estimated uncertainties (2-2.5%) 



tion (a! 



are represented by the band delineated by the dotted lines. 




0.37 



Figure 1: The pionium lifetime as a function of the combination 



of the S'-wave 



scattering lengths (full line). The band dehneated by the dotted lines takes into account 
the estimated uncertainties (2-2.5%). 



The theoretical value of the lifetime depends also on the uncertainties on the strong 
interaction threshold parameters, which were not considered above (these are not yet pub- 
lished in the literature in a definite form). A close analysis of the results obtained above 
shows, however, that eventual uncertainties on the threshold or low energy parameters 
have the tendancy to move the predicted value of the lifetime by remaining within the 
band of uncertainty already considered. An illustration of this phenomenon can be seen 
by considering the set II solution obtained in the standard X-P^ case |^ and for which 
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one has in particular (oq — a^) = 0.250. (A recent analysis of the threshold parameters 
in the standard case is presented in Ref. [0.) One might for instance consider 

the differences between predictions of set I and set II as uncertainties. One then has two 
possibilities of proceeding. First, one can calculate directly from set II the value of the 
lifetime, by repeating the calculations done with set I. Second, one can seek from the 
curve of Fig. |^ the value of the lifetime corresponding to (ctg — Oq ) = 0.250. For the first 
method, we have found from the threshold parameters of set II the following values of 
the low energy constants: a = 1.010, (3 = 1.111, Ai = -9.11 x 10"^, A2 = 11.15 x 10"^, 
A3 = -1.14 X 10"^, A4 = -0.27 X 10"^ The predicted value of the lifetime is r = 3.22 
fs. This value lies within the band of uncertainty of Fig. |^. The value one obtains for 
(oq — Oq) = 0.250 from Fig. |l], corresponding to the initial data of set I, is r = (3.20±0.07) 
fs, which underlines the fact that uncertainties of the order of 3% in the scattering length 
values have not moved the lifetime value outside the initial uncertainty band (and fur- 
thermore have left it rather close to the central line). Therefore, one should not enlarge 
further the initial uncertainty band of 2-2.5%, which could be considered as a conservative 
one. 

One thus arrives at the conclusion that the relationship between the pionium lifetime 
and the combination (oq — Oq) of the S'-wave scattering lengths represented in Fig. |1|, 
together with its uncertainty band, is practically independent of the particular choice of 
the low energy constants which yield the value of (cq — a^). In this sense, the curve of 
Fig. |T] should be used to determine from the experimental value of the pionium lifetime 
(with uncertainties) the corresponding value of (og — Oq) (with corresponding uncertain- 
ties). The determination of the other thereshold parameters or low energy constants, and 
in particular of a and of the quark condensate parameter, requires a separate analysis 
with the aid of additional constraints, coming for instance from sum rules and the Roy 
equations. 

It is also evident that for a resulting value of (oq — Oq), because of the other existing 
constraints, the value of a, and hence of the quark condensate parameter, could not be 
arbitrarily varied. Independent of any detailed analysis which should determine the precise 
value of a, one still is allowed to derive from Fig. |I| some qualitative conclusions which 
we summarize as follows. Values of the lifetime close to 3 fs, lying above 2.9 fs, say, would 
confirm the scheme of standard xPT. Values of the lifetime lying below 2.4 fs remain 
outside the domain of predictions of standard xPT and would necessitate an alternative 
scheme of chiral symmetry breaking. Values of the lifetime lying in the interval 2.4-2.9 
fs, because of the possibly existing uncertainties, would be more difficult to interpret and 
would require a more refined analysis; they might also indicate, within the standard xPT 
scheme, a strong dependence of the quark condensate on the light quark fiavor number, 
as a trace of a possible chiral phase transition at higher values of the latter. 
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